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Abstract 

Total and differential cross sections for the exclusive reaction pp ppp observed via the tt+tt" 
decay channel have been measured at Pbeam = 3.67 GeV/c. The observed total meson production 
cross section is determined to be (23.4 it 0.8 it 8)//6 and is significantly lower than typical cross 
sections used in model calculations for heavy ion collisions. The differential cross sections measured 
indicate a strong anisotropy (~ cos^O^^) in the meson production. 

PACS numbers: PACS numbers: 14.40.Cs, 13.75.Cs, 25.40 Ve 
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The study of the production of vector mesons (such as p^) in both hadronic and electro- 
magnetic processes is considered an excellent tool to investigate the properties of the hadrons 
both inside the nuclear medium and in free space. These properties are closely related to the 
QCD vacuum structure characterized by the presence of a non-vanishing quark condensate 
breaking chiral symmetry In QCD inspired models, in- medium p° mass modifications 
have been proposed as a signal for chiral symmetry restoration in dense and/or hot baryonic 
matter [jl| . Based on these calculations sizable changes in the p° spectral function have been 
predicted for near-threshold p° production in proton-, pion- and nucleus- nucleus reactions 

1^ . Furthermore modifications of the in- medium spectral functions can also rise from 
p — N coupling as indicated by hadronic model calculations @, 0. Such modifications could 
have a significant impact on the role of the as one of the mediators of the nuclear force 
at small distance, mainly in the tensor part of the interaction 0, |^. 

A measurement of the p production cross sections in pp collisions close to threshold is 
of particular importance for the understanding of the p — N coupling. Unfortunately, the 
large p width F = 0.15 GeV/c^ and a small cross section hampered up to now successful 
p° identification for excess energies above threshold e < 1.1 GeV In this letter, we 
report on the first measurement of the exclusive total and differential cross sections of the 
p° production in proton-proton reactions at e = 0.33 GeV, measured at the SATURNE II 
facility at Saclay. 

The proton beam of momentum 3.67 GeV/c was incident on a liquid hydrogen target 
and charged particles were detected using the DISTO spectrometer, which is described in 
detail elsewhere [|lOl- This spectrometer consisted of a large dipole magnet (40 cm gap 
size, operating at 1.0 T-m) which covered the target, located at the center, as well as two 
sets of scintillating fiber hodoscopes. Outside the magnetic field, two sets of multi-wire 
proportional chambers were mounted, along with segmented plastic scintillator hodoscopes 
(SH) and water Cerenkov detectors (CD) for particle identification. 

The acceptance of all detectors (~ 2° to ~ 48° horizontally and ~ ±15° vertically) on 
both sides of the beam, allowed coincident detection of four charged particles with sizable 
efficiency. The multi-particle trigger, which was based on the multiplicity of SH elements 
and the scintillating fibers, selected events with at least three charged tracks in the final 
state. With this trigger many exclusive channels, e.g. ppK~^K~ , ppTc^n^TT^ [ITT], ITEI and 
pKA, pKH |1^ were simultaneously measured. 
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The exclusive p° meson production {pp ppp^) was measured via its dominant (~ 100%) 
TT+TT^ decay channel. Proton and pion identification was achieved by means of conditions 
defined on two dimensional distributions of particle momentum versus Cerenkov light and 
energy loss measurements in the CD and the SH, respectively. It was found that less than 
10% of the pions could be misidentified protons. The selection of the pp pp-K^Ti^ reaction 
was based on two kinematical conditions that were checked with Monte Carlo simulations: 
(a) the four particle missing mass has to be zero {M^'^^Y = ± 0.015 GeV^/c^ and (b) 
the proton-proton missing mass has to be equal to the vr+vr" invariant mass (M™'^'')^ = 
{M^^+^-Y ± 0.2 GeV^/c"^. The residual background below the (M^*''*)^ bump amounts only 
to about 2% of the total yield and can mainly be attributed to events from the ppn^n^n'^ 
reaction. 

The acceptance of the DISTO spectrometer for the ppTr+vr" channel was determined using 
Monte Carlo simulations, which after digitization, were processed through the same analysis 
chain as the measured data. The production of a meson X with the given mass Mx, in the 
exclusive reaction pp —>■ ppX — > ppTT^n^ , was assumed for the final state particles generation. 
Eight linearly independent degrees of freedom are necessary to fully describe a four body 
final state. We have chosen the following variables: (M™)^)^, {M''^xY ^^'^ three Euler angles 
6x'^,(j)'x^,iljpj^^ in the center of mass frame (CM) that characterize the ppX system, and 
Mx,0^ ,(f)^ which describe the X —>■ tt+tt^ decay []T3|. The spectrometer acceptance has 
been determined as a four dimensional function of the variables {Mp^^Y, (M^"^)^, and 
Mx, assuming isotropic distributions in the remaining four variables. Symmetry reasons 
ensure that the distribution of (px^'^ must be isotropic, and the isotropic assumption for the 
remaining three variables was corroborated by the measured data. 

The simulations indicate a very low acceptance of the apparatus for the X produced in 
the backward hemisphere in the CM frame. However, since the initial system consists of 
two identical particles a refiection symmetry about 6^^'^ = 90° exists, thus the backward 
data are redundant for the total cross section determination. In the forward hemisphere 
the detector acceptance was found to be non-zero in every bin over the full kinematically 
allowed region, thus eliminating the need for any model dependent extrapolations of the cross 
section into unmeasured regions of phase space. As a result, the acceptance corrections for 
the DISTO spectrometer are independent of the actual final state distributions, except for 
residual effects related to the size of the kinematic bins and the detector resolution, which 
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FIG. 1: Proton-pion invariant mass distributions Mp™ selected according to the absolute value 
of the four-momentum transfer t: M*™+ with the lower | t \ value (left frame) and with the 

larger | t \ (right frame). The curves are fits to the data as explained in the text. 

have been accounted for in the estimate of the systematic error. Finally, the data presented 
below have been corrected on an event-by-event basis, via a weighting factor obtained from 
the acceptance function for the appropriate kinematic bin. 



Compared to the other mesons identified in our experiment [IT], |I2|, |l^] the extraction of 
the p° yield from the non resonant background is much more difficult because of its larger 
width. Therefore, different selection criteria had to be applied. For each of them the effect 
on the acceptance for and background events has been evaluated using the Monte Carlo 
simulation mentioned above. A first selection was applied on the three body pvr'^'vr" invariant 
mass Mp"+^_ > 1.6 GeV/c^, for both ptt+tt^ combinations. It was found that this rejects 
35% of the ppn^Ti^ events, whereas the simulations indicate that only less than 4% of p° 
events are lost by this restriction. 

A further investigation of the remaining events shows that, as observed at higher ener- 
gies [y], most two pion events originate from the pp A++p7r~ channel. Figure ^ shows the 



acceptance corrected pvT"*" invariant mass distributions for events fulfilling the M*"+^_ > 1.6 
GeV/c^ condition discussed above. The Mp™+ and Mp™+ have been obtained selecting the 
events according to the absolute value of the four momentum transfer \ t \, from the incom- 
ing beam proton to the outgoing p7r+ system. Mp™+ and M™^+ denote the invariant mass 
of the pvr^ pair with smaller and larger momentum transfer, respectively. The prominent 
A++ peak at Mj^^l = 1.231 ± 0.008 GeV/c^ is visible in both distributions. A fit to the 
data using a relativistic Breit-Wigner distribution P] of width F = 0.12 GeV/c^ added onto 
a third order polynomial background shape indicates that 77% of all ppir'^n' events involve 
^++ production. 

In order to reduce the background related to production, it was required that the 
invariant mass for both pn'^ combinations be 1.15 GeV/c^ < M™+ < 1.28 GeV/c^ and that 
the three body invariant mass Mp™+^_ < 1.74 GeV/c^. With such conditions 50% of the 
TT+TT^ pairs were rejected whereas the simulation showed that only 10% of those from p° 
decay are lost. A similar selection on pir^ invariant mass was judged unnecessary, since 
events involving A'^ production are much less abundant than those with A++ and therefore 
will not enhance the p° signal to background ratio. 

For the events surviving all the criteria described above, figure ^ presents the final, accep- 
tance corrected M^^^- distributions for four equal bins in cos {6'^^^) between cos(6 ) = 0.2 
and cos{6'"^^) = 1.0, where 6''"^^ is the CM polar emission angle of the di-pion system with 
respect to the beam axis. The dashed line shows the smoothed M^+'"^_ distribution for 
< cos{6^^^) < 0.2 normalized to the maximum, which reproduces well the general shape 
of the 7c'^'K~ non resonant spectra. It should be noted that this yield is almost independent 
of cos(^'"*^) indicating a nearly isotropic vr^vr^ non-resonant production. A structure near 
the p° mass Mp = 0.77 GeV/c^, is clearly visible. 

The p° yield has been extracted from a fit to the M^+"^_ distributions using a function 
a{m) = BW{m,r) ■ V^s(m) describing the line shape and a third order polynomial 
background. The BW{m, F) is a relativistic Breit-Wigner function with a mass dependent 
pion decay width of the meson given by P]: F = F(mij) • (y^^^^)'^ • where F(mij) = 0.15 
GeV/c^ is the full p° width at its pole position m^j = Mp and k.,^{m), A;^(mR) are pion three- 
momenta in the rest frame of the resonance with mass m and m^, respectively. Vps{m) 
stands for the phase space available for the direct (no intermediate state) production of the 
meson with a given mass m. We have calculated Vps{m) by numerical integration of the 



6 



150 



100 



50 



150 



100 



50 



1 1 1 1 1 1 1 1 
ii'ii 'A i! 

/' \ 

O.2<cos(0'")<O.4 \ 

1 1 


1 1 1 1 1 1 1 1 1 

> 

O.4<cos(0''')<O.6 \ 

1 1 1 1 1 1 1 1 1 ' 


/ v., 

• V 

O.6<cos(0'")<O.8 \ 


/ 

■ » .* 

h \ 

'I 

O.8<cos(0™)<1 \ 



0.4 0.6 OA 



0.4 



0.6 0. 

inv 

'tttt 



[GeV/cl 



FIG. 2: vr^vr^ invariant mass distributions obtained in four different CM angular bins as indicated 
on the figure. The dashed hne represents the smoothed invariant mass distribution for the < 
cos{6^^^) < 0.2 normahzed to the maximum. 

simulated ppp Dalitz distribution including the suppression cuts discussed above. As 
an example, figure |^ shows the M*+'^_ distribution for cos{6'-^^'^) > 0.5 containing most of 
the p° yield. The initial background parameters have been obtained by fitting the M^^- 
distribution for cos(6''"^^) < 0.4 and M^™^_ > 0.55 GeV/c^ normalized to the same maximum 
(dashed line). The best fit ij^ jv = 1.2) has been obtained with a fixed Tlnifj) = 0.15 and 
Mr = 0.772 ± 0.003 GeV/c^ (dotted) and background (solid line) which differs only very 
little from its initial shape. 
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In order to extract the angular distribution for the production and to extract the total 
production yield, a finer angular bin size was chosen and the fit procedure outlined above 
has been applied. The total p° yield has been obtained as a sum of yields in the individual 
angular bins. The total cross section of a = (23.4 ± 0.8 ±8)yu6 with statistical and systematic 
errors, respectively, was determined by normalizing the acceptance corrected p'' yield to the 
simultaneously measured exclusive r] yield of known cross section o"^ = 135 ± 35 /ib. The 
T] yield was obtained using separate acceptance correction matrices but calculated with the 



same method |13|. This procedure to determine the absolute normalization was chosen 
because it reduced the large systematic uncertainty associated with the determination of 
the absolute beam intensity and trigger efficiency. It was checked by detailed Monte Carlo 
simulations that particle identification efficiencies and yield reductions due to the analysis 
cuts applied for the selections of various reactions are understood within 15% [0. Finally, 
the method was verified by comparing the exclusive cross section for uo production measured 



in our experiment [|T3| with other experimental results obtained in a similar energy range. 
An interpolation between existing data points give 0"^^ = 45 ± 7 p6 |jl6| in good agreement 
with our result (To; = 50 ± 3lj;gp6. 

The systematic uncertainty for the p° cross section contains two almost equal contribu- 
tions: (a) errors from absolute normalization and (b) uncertainties related to the background 
subtraction, acceptance corrections, tracking efficiency and analysis cuts. In particular a pos- 
sible bias in the analysis due to the A^"*" suppression has been investigated by comparing 
the p° intensity obtained with and without the suppression. The observed yield reduction 
of ~ 10% is reproduced by the Monte Carlo simulations mentioned above, was corrected for 
and completely included in the systematic errors. 

Polar angle differential cross sections have been proposed as a sensitive observable to 
extract the strength of nucleonic and mesonic currents responsible in meson exchange models 
for the vector meson production near the threshold. Therefore, it is particularly interesting 
for various reaction models to compare the differential cross sections for the production 
of p'' and uo as measured in our experiment |13[. The lower part of Figure |^ shows the 
p*^ (circles) and uj (triangles) production cross section as a function of cosO'"^^ , the polar 
angle in the CM system of the meson. The vertical error bars reflect the statistical errors, 
only. The distributions have been flt with the sum of the flrst three Legendre Polynomials 
(solid line) Pi with the result: ^ oc (2. ± 0.2)Po + (4.4 ± 0.2)P2 + (0.8 ± 0.2)P4 and 
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^ oc (4 ± 0.1)Po + (3.1 ± 0.2)P2 + (2.0 ± 0.2)P4 for the p° and respectively. 

The above results indicate contribution of partial waves up to L=2 in the production 
of both mesons. Clearly, the angular distribution for p° production exhibits stronger for- 
ward peaking, suggesting a dominance of the nucleonic current characterized by a dipole 
~ cos^O^^'^ dependence The angular distributions for uj production at this beam energy 



and at the smaller excess energy e = 0.173 GeV |[T^ show strong nucleonic contribution too, 
but also an additional isotropic component signaling the importance of the mesonic np ^ u 
fusion It is also interesting to note, that the observed ratio of the total cross sections 
for p° and uj production R = (y{p^) / a{uj) = 0.46 ± 0.12 is very similar to the ratio obtained 
at higher energy (e = 1.14 GeV) R = 0.48 ± 0.13 || and seems to indicate similar scaling 
of the cross section with -^i for both mesons. 

In conclusion, the production of the p° meson has been studied in the pp reaction at 
Pbeam = 3.67 GcV/c. A total p° cross section of cr = (23.4 ± 0.8 ± 8)p6 has been determined. 
It is significantly lower than typical cross sections used by the hadronic models (~ 100 fib) 
for meson production in pA and AA collisions 0, [l^. Furthermore, the differential cross 
sections exhibit a prominent ^ oc dependence that might signal the dominance of 

nucleonic current in p^,^^ meson production. 

We thank the SATURNE II staff for delivering an excellent proton beam and for assisting 
this experimental program. 

This work has been supported by: CNRS-IN2P3, CEA-DSM, NSF, INFN, KBN 
(5P03B14020) and GSI. 

Now at: DAPNIA/SPhN, CEA Saclay, France. 
LPHNHE, Ecole Polyt. 91128 Palaiseau, France. 
^ Temple Univ., Philadelphia, Pennsylvania, USA 

Paul Scherrer Institut, Villigen, CH-5232. 
^ Motorola Polska Software Center, Krakow, Poland. 
^ Dip. di Fisica "A. Avogadro" and INFN- Torino, Italy. 
^ Arthur Andersen, Eschborn, Germany. 
^ lU School of Medicine, Indianopolis,USA 



9 



[1] G. E. Brown and M. Rho, Phys. Rev. Lett. 66, 2720 (1991). 

[2] G. Q. Li,C. M, Ko, and G. E. Brown, Phys. Rev. Lett. 75, 4007(1995). 

[3] W. Gassing and E. Bratkovskaya, Phys. Rep 308, 65 (1999). 

[4] R. Rapp and J. Wambach, Adv. Nucl. Phys. 25, 1 (2000). 

[5] W. Peters et al., Nucl. Phys. A632, 109 (1998). 

B. Friman and H. Pirner, Nucl. Phys. A617, 496 (1997). 

[6] R. Machleidt,Advances in Nuclear Physics 19, 189 (1989) 

[7] D. Jido, E. Oset and J.E. Palomar, |nucl-th/0202070| 

[8] K. Nakayama et al., Phys. Rev. G 57, 1580 (1998). 

[9] G. Alexander et al., Phys. Rev. 154, 1284 (1967). 
E. Golton et al., Phys. Rev. D 3, 1063 (1971). 
[10] F. Balestra et al., Nucl. Instr. Meth. A 426, 385 (1999). 
[11] F. Balestra et al., Phys. Rev. Lett. 81, 4572 (1998). 
[12] F. Balestra et al., Phys. Lett. B468, 7 (1999). 
[13] F. Balestra et al., Phys. Rev. G 63, 024004 (2001). 
[14] F. Balestra et al., Phys. Rev. Lett. 83, 1534 (1999). 
[15] S. Teis et al., Z. Phys. A 356,421 (1997). 

G. Ernst et al., Phys, Rev. G 58, 447 (1998). 
E. Bratkovskaya, Nucl. Phys. A 653, 301 (1999). 
[16] A. Sibirtsev and W. Gassing, Eur. Phys. J.A7,407(2000). 
[17] S. Abd El-Samad et al., Phys. Lett B 522, 16 (2001). 



10 




0.4 0.5 0.6 0.7 0.8 0.9 

M':: [GeV/c'] 




0.2 0.4 0.6 0.8 1 

cos(0'=") 



FIG. 3: Upper frame: M^+j^- distribution measured for cos{9 ) > 0.5. The solid line shows the fit 
using a relativistic Breit-Wigner function added onto a third order polynomial for the background, 
as described in the text. The dotted line showg^|he meson line shape and the dashed line presents 
the M^"^_ distribution obtained for cos{9'^'^) < 0.4 normalized to the same maximum. Lower 



